Europhysics Letters 



PREPRINT 



Structural distortions and orbital ordering in LaTiOs and 
YTiOs 

S. Okatov^^^, a. Poteryaev^'' and A. Lichtenstein^ 

^ Research Institute for Materials, University of Nijmegen, The Netherlands 

^ Institute of Solid State Chemistry, UrD of RAS, Ekaterinburg, Russia 

^ I. Institute of Theoretical Physics, University of Hamburg, Germany 

^ Centre de Physique Theoretique, Ecole Poly technique, Palaiseau, France 



PACS. 71 .20.-b - Electron density of states and band structure of crystalline solids. 
PACS. 75.25.+Z - Spin arrangement in magnetically ordered materials. 
PACS. 71.27.+a - Strongly correlated electron systems; heavy fermions. 



Abstract. - Theoretical investigations of the electronic, magnetic and structural properties 
of LaTiOs and YTiOa have been made. In the framework of GGA and GGA+U scheme we 
analyzed the effect of the local Coulomb interaction [U) value on the atomic forces acting in the 
experimental structure. The optimal parameters of the electron-electron on-site interactions as 
well as the orbital configurations and magnetic properties are determined. 



Introduction. - Transition-metal perovskite oxides LaTiOs and YTiOa are classical Mott- 
Hubbard insulators. In spite of the fact that they are formally isoelectronic with one 3d 
electron in the t2g shell LaTiOa has a G-type antiferromagnetic (AFM) structure whereas 
YTiOa has a ferromagnetic (FM) one. The most unusual feature of these compounds relates 
with a nearly isotropic magnon spectrum in both titanates despite of their distorted crystal 
structures [1,2]. 

There are two theoretical models aimed to explain these effects in LaTiOs. One of them 
originates from the work of Kugel and Khomskii [3] and operates in terms of the lattice 
distortions and an orbital ordering (00). Another one is proposed by Khalliulin and Maekawa 
[4] and describes the above effects by the model of orbital liquid (OL) where the fluctuations of 
orbital degrees of freedom play an important role. According to [4] the energy of the orbital 
fluctuations is estimated to be Worb ^ 160 meV. Thus OL model is applicable only if the 
splitting of the t2g orbital in the crystal field A is significantly smaller than Worb- Otherwise 
the nearly isotropic exchange is likely to be due to an orbital ordering with a peculiar orbital 
configuration [9]. 

The OL model was supported by neutron experiments of Keimer et al [1] where no 00 
in LaTiOa was found. Fritsch with co-workers [5] made heat capacity and magnetic measure- 
ments. Having supposed the nearly cubic structure of LaTiOa they concluded that the energy 
of the spin-orbit coupling for Ti-sites was E^q = 30 meV and the crystal field effect was small. 
The theoretical estimations made by Solovyev in Ref. [6] in the framework of the local spin 
density approximation (LSDA) and LSDA-I-U theory showed that A in LaTiOa was as small 
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as 49 meV and this value was of the same order of magnitude as the energy of the spin-orbit 
(SO) coupUng £;so = 23 meV. 

Recent X-ray studies [7, 8] revealed a sizable deformation of the TiOe octahedra and de- 
viation from the cubic symmetry in LaTiOa below Tn. The splitting of the t2g- levels due to 
crystal deformations was estimated in [7] to be about 240 meV, and 120-300 meV in [9]. This 
was also supported in the ab-initio calculations [10] in the framework of LDA and LDA-I-DMFT 
theory, where A was reported to be ^ 200 meV. Those values of A were larger than Worb and 
could provide ordering of the t2g orbitals. The 00-model has had a support in the ab-initio 
theoretical studies [10-13] for both YTiOa and LaTiOs as well as in the direct observations 
of the orbital ordering in YTiOs [14-16], and in LaTiOs [16] from the NMR spectra analysis. 

The physical properties of both LaTiOs and YTiOs are strongly sensitive to the chemical 
composition [17]. The highest value of the Neel temperature in LaTiOa was reported by Cwik 
et al [7] and Hemberger with co-authors [8] to be Tn = 146 K, while the magnetic moment was 
IJ, = 0.57/iB [7] . This magnetic moment disagreed with the theoretical estimations of more than 
0.8/ZB [6, 18]. For YTiOs the Curie temperature amounts to 30 K and fj, = 0.84/xb [15, 16, 19]. 

According to the 00 model the difference in the magnetic properties of LaTiOs and 
YTiOa can be attributed to the different types of distortions (see for example Refs. [6,7,13]). 
The main distortion in YTiOa has the Jahn- Teller type, while the GdFeOa-type distortion 
dominates in LaTiOa. In this letter we investigate the nature of those crystal distortions and 
the effects of 00 on electronic and magnetic properties of the LaTiOa and YTiOa taking the 
effects of local Coulomb correlations into account. 

Methods and models. - Both LaTiOa and YTiOs have a Pbnm crystal structure. The 
lattice constants and atomic positions are taken from [7,19]. The structural parameters for 
LaTiOa used in the calculations correspond to the temperature below Tn (T = 8 K). The 
examination of the YTiOs is made for the high temperature structure (T — 293 K) only 
because there is no structural data for YTiOs at T < Tc. The crystal structure relaxations 
are not performed but the calculated forces on atoms for the experimental atomic positions are 
used for the verification of the agreement between theory and experiment. The orthorhombic 
unit cell used in our calculations consists of four formula units. This cell allows us to consider 
both FM and AFM structures with different types of the magnetic ordering. 

For electronic structure calculations we use the projected augmented wave (PAW) method 
[20,21] in the framework of the density functional theory (DFT). The exchange-correlation 
correction is taken into account within the general gradient approximation (GGA) [22]. In 
order to treat the effect of local Coulomb interactions in a partially filled 3d band of Ti the 
LDA-l-U method is apphed [23,24]. The Brillouin zone (BZ) is sampled with the 7 x 7 x 5 
mesh with its origin at F point. The valence states include 2s, 2p for O, 4p, Ad, 5s for Y, 3p, 
3d, 4s for Ti and bp, bd, 6s for the La-atoms. The cutoff energy of the plane-wave expansion 
is 400 eV. 

Results and discussions. - The results of our GGA calculations for LaTiOs and YTiOa 
are presented in fig. ^ We have found the non-magnetic ground state for LaTiOs and fer- 
romagnetic ground state for YTiOa with the value of the magnetic moment O.S/ie per Ti 
atom. This magnitude of the magnetic moment is very close to the experimental value of 
0.85/iB [15,16,19]. 

In the right part of fig. ^ the total density of states (DOS) for ferromagnetic YTiOa 
calculated within GGA approach is shown by red color. It consists of three well separated 
sets of bands. The lowest in energy band lies from -8 till -4 eV and is predominantly of 0-2p 
character. The t2g states of Ti are strongly polarized. Spin up crosses the Fermi level whereas 
the spin down is completely empty and lies above the Fermi energy, thus YTiOa is a half 
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magnetic ferromagnet in our GGA calculations. The bandwidth of the t2g states is about 
1.8 eV for both spin channels. The eg states of Ti and the rest unoccupied states (03s, Ti4p, 
etc.) are higher in energy and separated from t2g states by a small gap. Our results are in 
good agreement with the early study of this compound [12,25]. 

In the upper left part of fig. ^ the total DOS for LaTiOa is shown. In this case the DOS 
is non-spin-polarized and consists of four parts. As in YTiOa the fully occupied band of 0-2p 
character is located from -7.8 to -4 eV. The t2g states cross the Fermi level and have the 
bandwidth about 1.9 eV. The eg states lie starting from 2 eV. In the energy region between 
t2g and eg states the La-4/ band is located. This band ordering Ti-t2g, La-4/, Ti-Cg is not 
compatible with the results of early calculations [6, 10, 11, 18,26] and with the experimental 
data [27]. The simplest way to fix this problem is to apply GGA-I-U approach for the La-4/ 
states. The results of the calculation of LaTiOa with a value of Coulomb repulsion U^'j = 8 eV 
are shown on the left down part of fig. Q by blue color. The value of Coulomb interaction 
used is in good agreement with the experimental and theoretical estimations [11, 18,28]. The 
obtained band ordering is now in consistency with the previous studies [10]. In the rest part 
of the paper we will always use this result as a starting point for the following considerations 
without explicit mentioning and will refer to this calculation as GGA. 

According to experiments [27,28] both LaTiOa and YTiOs are classical Mott-insulators 
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with the band-gaps (AE'g) of 0.2 and 1.2 eV respectively whereas in our GGA calculations 
they are metals. We have found also that the forces acting on the oxygen atoms in the GGA 
are quite big 0.5 eV/A) for the experimental crystal structure. 

It is well known that in the systems where the value of kinetic energy is comparable 
with the value of one site Coulomb interaction and strong electron-electron correlations take 
place, GGA (or LDA) approximation may fail in the prediction of the correct ground state. 
The reason is an underestimation of the strong correlations in DFT. In order to take into 
account the one site Coulomb interaction we use the LDA+U approach which has orbitally 
^m-dependent potential and describes well the insulating compounds with the long range 
magnetic ordering. In the LDA-I-U approach the potential is defined [23,24] as 

^LSDA+U ^ ^LSDA ^ _ ^) [ 1 ^^.^ _ j 

where pjj is the density matrix of d electrons. 

The value of the Coulomb parameter Uj^ was varied from 1 up to 5 eV. Since the value of 
the intra-atomic Hunds rule exchange parameter is almost independent on materials [23] we 
use a fixed value J — 1 eV. The results of the GGA-I-U calculations for LaTiOa and YTiOa 
with eV are shown by black color in fig.^ In both cases the oxygen bands are slightly 

higher in energy ~ 0.1 eV in comparison with the GGA results. The t2g band is splitted and 
the value of the gap is equal to 0.4 eV and 0.6 eV for LaTiOa and YTiOs respectively. The 
magnetic moment per Ti ion in LaTiOa grows fast with the value of Coulomb interaction 
and saturates to 0.83/iB at ~ 4 eV (see fig. The experimental magnitude of the 
magnetic moment for LaTiOa corresponds to Uj^ ~ 2 eV in our calculations. The insulating 
gap appears at Uj^ ~ 3 eV. In contrast to the GGA calculation the ground state of LaTiOs 
is the G-type AFM with the energy gain relative to other types of magnetic ordering of 68 
(FM), 85 (C-type AFM) and 23 meV (A-type AFM). In YTiOs the value of gap increases 
with the magnitude of the Coulomb interaction (not shown). The experimental gap of 1.2 eV 
wide corresponds to parameter Uj^ — 4. The magnetic moment in YTiOa (fig. ^ is almost 
constant and equal to ~ 0.8 eV. The ferromagnetic solution is the lowest in energy. 

In fig. [2 we plot the modules of the forces acting on 01 and 02 atoms versus the value of 
the Coulomb interaction. One can see the forces for 01 and 02 in LaTiOs and 02 in YTiOa 
reduce with increasing of U^^ and reach their minimum at about 3.5 cV. The forces for 01 in 
YTiOs increase with Ujj. This fact can be related to the high temperature structure used, 
thus we expect large distortions in YTiOa below Tq. The value of Uj^ = 3 — 4 eV is in good 
agreement with the previous experimental [2,28] and theoretical [18] estimations. 

The directions of the forces calculated in the GGA and GGA-I-U {Uji^ — 3 eV) approx- 
imations are shown in fig. 01 by yellow and blue colors respectively. In the case of the GGA 
method one can see that in LaTiOs they are almost equal in absolute value and trend more 
to rotate the octahedra than to distort them. In YTiOs the two inward forces face nearly 
along Ti-0 bonds, while the two rest outward ones rotate. Thus the GGA method tries to 
make the structure more cubic. That is a characteristic of the titanates at temperatures above 
the orbital ordering (700 K [8]) where the electrons become more delocalized and are better 
described by the GGA. In the case of GGA-f U the forces smoothly decrease and rotate with 
increasing of Uj^ (not shown). At Uj^ = 3 eV the angle between GGA and GGA-I-U forces is 
- 90° in LaTiOa and ^ 30° in YTiOg (fig. 13 . 

In order to visualize the resulting orbital ordering we plot the charge density in fig. |2| for 
both titanates integrated in the energy window from —2 eV till Fermi level. In this energy 
interval the charge density has predominant Ti-t2g character. One can see different types of 
ordering in LaTiOa and YTiOa: whereas in LaTiOa the orbitals arrange in a "fish bone" style 
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they form a rhombus in YTiOs. These results agree with those reported in previous theoretical 
findings [10] and experimental works [7,14-16]. The presented electronic densities clarify the 
mechanism of the lattice distortions. Because of the repulsion between negatively charged 
orbitals and oxygen atoms the different types of 00 cause diffcrnt types of the distortions. 
The orbitals in LaTiOa do not face toword any oxigens, thus distortions of the octahedra are 
small there and their rotations are dominant. On contrary in YTiOa a specific elongation of 
the orbitals toword 0-atoms makes strong distortions of the octahedra favourable. Since 00 
in the GGA calculations is weak the corresponding forces face face opposite direction. Thus 
taking into account the one site electronic correlations with the energy of Coulomb interaction 
around 3 — 4 eV we obtain a correct type of the orbital ordering and can describe the crystal 
distortions. 



LaTiOs YTiOs 




Fig. 3 - {color online) Charge densities in the energy window from -2 eV till EF (gray surface). 
Arrows are forces acting on atoms in GGA (orange) and GGA+U with Ujj^ = 3 (blue). = 8 eV. 
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Conclusions. - We have found the correct magnetic and electronic ground state for 
LaTiOs and YTiOa in the framework of the GGA+U calculations. The value of Coulomb 
interaction parameter estimated from the minimum of the forces is about 3 — 4 eV. This is 
in good agreement with the previous theoretical studies and photoemission experiments [11]. 
In YTiOs in the GGA+U method with the optimal i7-values the band-gap is 0.6 — 1.2 eV 
and the magnetic moment is ^ = 0.81 — 0.83/xb in a good agreement with experimental data. 
The values of /x = 0.7/zb and gap Eg = 0.4 eV for U^J^ = 3 eV in LaTiOs are overestimated 
because of the absence of the magnetic fluctuations in GGA+U [10]. The orbital order found 
for all values of Uj^ is "fish-bone" and rhombus-type for LaTiOs and YTiOs respectively and 
is consistent with experimental structure findings. 
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